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Abstract 

We consider the effect of Planck scale operators on neutrino mixing. We assume that GUT scale 
operators give rise to degenerate neutrino masses with bimaximal mixing. Quantum gravity (Planck 
scale) effects lead to an effective SU (2)^ x [/ (1) invariant dimension-5 Lagrangian involving neutrino 
and Higgs fields. This gives rise to additional terms in the neutrino mass matrix on electroweak 
symmetry breaking. These additional terms can be considered as a perturbation to the GUT scale 
bi-maximal neutrino mass matrix. We assume that the gravitational interaction is flavour blind 
and compute the deviations of the three neutrino mixing angles due to the Planck scale effects. We 
find that the changes in ^13 and 623 are very small but the change in solar mixing angle 9i2 can be 
as large as 3.5°. 
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I. INTRODUCTION 



One of the challenges in neutrino physics is to explain the pattern of neutrino masses and 
mixings which are deduced from the current neutrino data. Presently accepted values of 
the mixing angles are 612 = 34" ± 7°, ^13 < 12° and 623 = 45° ± 9° It is reasonable to 
suppose that the symmetries of the Lagrangian predict the elements of MNS matrix to be 
simple fractions. In a well-known scenario called bi-maximal mixing j2|, the mixing angles 
are predicted to be 6^12 = 45°, 613 = 0° and 6^23 = 45°. Additional effects can modify the 
above predictions and lead to values of mixing angles close to experimentally determined 
values. For example, it is well-known that renormalization group evolution from GUT scale 

n 

to electroweak scale can substantially change the value of 6*12 [3||. The neutrino mass matrix 
is assumed to be generated by the see saw mechanism j^. Here we will assume that the 
dominant part of neutrino mass matrix arises due to GUT scale operators and they lead to 
bi-maximal mixing. This matrix can receive corrections due to physics from higher scale. 

The gravitational interaction of neutrinos with the Standard Model Higgs field can be ex- 
pressed as an effective SU{2)l x U{1) invariant dimension-5 operator [5|, 

Lgrav = ^ilpAaeAci^c)C~fl-ilpBbf}(^BDi'D) + h.C. (1) 

Here and everywhere below we use Greek indices a,/?. .for the flavour states and Latin 
indices i, j, k for the mass states. In the above equation ipa = {^a, la) is the lepton doublet , 
= 0°) is the Higgs doublet and Mpi = 1.2 x lO^^GeV is the Planck mass. A is a 3 x 3 
matrix in flavour space with each element 0(1). In eq(l), all indices are explicitly shown. 
The Lorentz indices a, b = 1, 2, 3, 4 are contracted with the charge conjugation matrix C and 
the SU(2)L isospin indices A, B, C, D = 1, 2 are contracted with e, the Levi-Civita symbol 
in two dimensions. After spontaneous electroweak symmetry breaking the Lagrangian in 
eq(l) generates additional terms to the neutrino mass matrix 

-1 

-'-'mass r 

-XapVaC Up, (2) 

where v=174 GeV is the VEV of electroweak symmetry breaking. 

We assume that gravitational interaction is "flavour blind", that is A^/? is independent of 
a, j3 indices. This is a reasonable assumption because gauge invariance requires gravity 
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to couple to the spin of a particle and these couplings are independent of the global U{1) 
charges of the particle. Thus the Planck scale contribution to the neutrino mass matrix is 

/l 1 l\ 



II X — /I 



111 
111 



(3) 



where the scale /j, is 



2.5 X IQ-^eV. 



pi 



(4) 



In our calculation, we take eq(3) as a perturbation to the main part of the neutrino mass 
matrix that is generated by GUT dynamics. We compute the changes in neutrino mass 
eigenvalues and mixing angles induced by this perturbation. 

II. CORRECTIONS TO MIXING ANGLES AND NEUTRINO MASS SQUARED 
DIFFERENCES 



We assume that GUT scale operators give rise to the light neutrino mass matrix, which in 
the mass eigenbasis, takes the form M = diag{ Mi,M2, M3), where Mj are real and non 
negative. We take these to be the unperturbed (0*'* — order) masses. Let U be the neutrino 
mixing matrix at 0*'* — order. Then the corresponding 0*'* — order mass matrix M in flavour 
space is given by 

M = U*MUl (5) 



The 0*'* — order MNS matrix U has the form 5 



I 



(6) 



^el ^e2 ^e3 
t^/il t^/i2 t^/iS 
V C/ri C/t-2 C/r3 } 

where the nine elements are, in general, functions of three mixing angles and six phases. In 
terms of the above elements, the mixing angles are defined by 

(7) 

23, (8) 
13- (9) 



C/e2| 


— tanO 


Uel 






— tanO 






\UA 


— sinO 



In terms of the above mixing angles, the mixing matrix is written as 

U = diag{e'^\ e'^\ e'^^)R{e23)AR{en)A*R{ei2)diag{e'''\ e'^^ 1). (10) 

The matrix A = diag{e~ , l,e~) contains the Dirac phase 6. This phase leads to CP vio- 
lation in neutrino oscillations, al and a2 are the so called Majorana phases, which affect 
the neutrinoless double beta decay. /I, /2 and /3 are usually absorbed as a part of the 
definition of the charged lepton field. It is possible to rotate these phases away, if the mass 
matrix in eq(5) is the complete mass matrix. However, since we are going to add another 
contribution to this mass matrix, these phases of the zeroth order mass matrix can have an 
impact on the complete mass matrix and thus must be retained. By the same token, the 
Majorana phases which are usually redundant for oscillations have a dynamical role to play 
now. Planck scale effects will add other contributions to the mass matrix. Including the 
Planck scale mass terms, the mass matrix in flavour space is modified as 

M ^ M' = M + /iA, (11) 

with A being a matrix whose elements are all 1 as discussed in eq(3). Since fi is small, 
we treat the second term (arising from the Planck scale effects) in the above equation as a 
perturbation to the first term (the GUT scale mass terms). This perturbation formalism 



was first developed in ref. 



6||. Here we briefly recall the main features for completeness. 



The matrix relevant for oscillation physics is the following hermitian matrix 

M'^M' =(M + fiX)\M + fiX). (12) 

To the first order in the small parameter /i, the above matrix is 

M^M + /iA^M + M^A. (13) 

This hermitian matrix is diagonalized by a new unitary matrix U' . The corresponding di- 
agonal matrix M'\ correct to first order in /i, is related to the above matrix by U' M''^U'\ 
Rewriting M in the above expression in terms of the diagonal matrix M we get 

U'm'^U'^ = U{M^ + m)M + Mm)U^ (14) 

where 

m = nU'XU. (15) 



Here M and M' are the diagonal matrices with neutrino masses correct to 0*'' and l*'^ order 
in /i. It is clear from eq(14) that the new mixing matrix can be written as: 

U' = U{l + iSe), (16) 

where 6Q is a hermitian matrix that is first order in /i. Oscillation physics is unchanged 
under the transformation U ^ UP, where P is a diagonal phase matrix. We can use this 
invariance to set the diagonal elements of the matrix 6Q to be zero. 

From eq(14) we obtain 

+ m^M + Mm = M''' + [i5Q,M"]. (17) 
Therefore to first order in /i, the mass squared difference AMfj = Mf — Mj get modified as 

AM''j = AMfj + 2{MiRe[mu] - MjRe[mjj]). (18) 
The change in the elements of the mixing matrix, which we parametrized by 6Q, is given by 

~ AMI AMI ■ ^^^^ 

The above equation determines only the off diagonal elements of matrix SQij. As mentioned 



ity. The expressions 
From the above 



above, the diagonal elements can be set to zero without loss of genera 
for AM -J in eq(18) and for 6Qij in eq(19) were first obtained in ref 
equation, we see that 6Qij are proportional to the neutrino masses. Thus they are larger 
for the case of degenerate neutrinos. Here onwards we assume degenerate neutrino masses. 
Then expression for 6Qij simplifies to 

In obtaining the above equation, we made the approximation AM^^ = AMf^. This is valid 
because the rriij term in the numerator is already proportional to /i and we are working to 
1** order in /i. For degenerate neutrinos, Mj + Mj ^ M^ — Mj. Thus the expression for SQtj 
further simplifies to 

^0.- = (21) 



To obtain the largest effect possible, we assume the largest allowed value of 2 eV for degener- 
ate neutrino mass which comes from tritium beta decay {t]. We also assume normal neutrino 



mass hierarchy. Thus we have Mi =2 eV, M2 = + A21 and M3 = ^Mf + A31. If all 
elements of A consist of 1 then rriij — ZiZj, where three complex numbers Zi are defined by 

Zl = Uel + U^l + UrU (22) 

Z2 = Ue2 + U^2 + Ur2: (23) 
Z3 = UeS + U^S + Ur3. (24) 

As in the case of O*'*order mixing angles, we can define 1''* order mixing angles in terms of 
1** order mixing matrix elements in a manner similar to eq(19) 

U' 

tan0[2=\jf'\ (25) 



el 



From eq(16), we get 



tan4 = &| (26) 
sin9[^^\U'j. (27) 



SUaj = Kj - U^j = i E Uai^Bi^ (28) 

i=l 



Substituting the expressions for 8Qij from eq(21) (and dQa = 0) in the above equation, we 
obtain 

5U,i = /. ( U.2^^^ + U.s^^^] , (29) 

For degenerate neutrinos, M3 - Mi = M3 - M2 > M2 - Mi because A31 = A32 > A21. 
Thus, from the above set of equations, we see that 5Uei and SUe2 are much larger than 
SUe3, SUfj,3 and 6Ur3. Hence we can expect much larger change in ^12 compared to ^13 and 

^23- 
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III. NUMERICAL RESULTS 



As mentioned in the introduction, we expect the mixing angles coming from GUT scale 
operators to be determined by some symmetries. We assume these symmetries give riso to 
a bi-maximal mixing pattern, 612 = 623 = vr/4 and 613 = 0. Using eq(25) to eq(27), we 
compute the modified mixing angles for the degenerate neutrino mass of 2 eV. We have 
taken A31 = 0.002eU2 (3| and A21 = O.OOOOSeU^ 9|. For simplicity, we have set the charged 
lepton phases /I = /2 = /3 = 0. Since we have set 613 = 0, the Dirac phase 6 drops out 
of the O*'^ order mixing matrix. From eq(22) to eq(24) we see that the complex numbers 
zi, Z2 and Z3 are independent of b and hence all deviations in U^j are independent of 5. 
The complex nature of Zi, Z2 and z^ comes from the Majorana phases al and a2, which we 
take to be non-zero. We compute the modified mixing angles as function of al and a!2. In 
table 1, we list the modified neutrino mixing angles for some sample values of al and a2. 
As shown in the table, the deviation in ^13 and 6*23 is negligibly small whereas the deviation 
in 6^12 is significant. In fig(l), we show contours of constant deviation, bQ\2 = 9'i2 ~ ^12), vs 
al and a2. 
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Figure 1: Deviation of 612 as a function of the Majorana phases in the case of bi-maximal mixing. 
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Figure 2: Deviation of Ou as a function of the Majorana phases in the case of tri-bi-maximal mixing. 



From fig(l), we see that Planck scale effects reduce 9i2 from the bi-maximal value of 45° to 
about 6*^2 = 41.5°. This is the present 3a upper limit of the solar mixing angle. Thus we see 
that Planck scale effects can bring down bi-maximal 612 to within experimentally acceptable 
range. In our discussion above, we have not included the renormalization group running. 
The matrix in eq(3) arises at Planck scale. Its renormalization group running from Mpi to 
Mgut is negligible. Thus, at Mqut, the total interaction matrix, that eventually gives rise 
to neutrino mass matrix on electroweak symmetry breaking, is the sum of GUT scale matrix 
and the Planck scale matrix. Renormalization group running modfies these couplings from 
their GUT values to low energy values ^. This can lead to lower values of 6'i2, much closer 
to the experimental best fit value. 



Recently various authors considered tri-bi-maximal mixing scenario which predicts sin'^Ou = 
1/3 or ^^12 = 35°, in order to obtain a value close to the best fit value lOj. We computed 
the modified mixing angles for this scenario also, that is for the input values 6*12 = 35°, 6*13 = 
0°and 023 = 45°. Neutrino masses are kept to be the same as in the previous case. The 
modified mixing angles are shown in table (2). In fig (2) the deviation of 6'i2is plotted as a 
function of al and a2. We see that deviation in 612 is about ±3°. The range of 6^12 given by 
this deviation spans la range of solar neutrino mixing angle. 
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al 


a2 


^12 


^23 


^13 








48.57 


45.00 


0.28 





45 


46.87 


44.93 


0.22 





90 


44.99 


44.86 


0.13 





135 


46.94 


44.93 


0.22 





180 


48.57 


45.00 


0.28 


45 





46.68 


45.07 


0.22 


45 


45 


44.96 


45.00 


0.20 


45 


90 


43.09 


44.93 


0.09 


45 


135 


45.03 


45.00 


0.14 


45 


180 


46.68 


45.07 


0.22 


90 





45.00 


45.14 


014 


90 


45 


43.28 


45.07 


0.10 


90 


90 


41.42 


45.00 


0.00005 


90 


135 


43.00 


45.07 


0.10 


90 


180 


45.00 


45.14 


0.14 


135 





46.68 


45.07 


0.22 


135 


45 


44.96 


45.00 


0.14 


135 


90 


43.09 


44.93 


0.09 


135 


135 


45.03 


45.00 


0.20 


135 


180 


46.68 


45.07 


0.22 


180 





48.57 


45.00 


0.22 


180 


45 


46.87 


44.93 


0.22 


180 


90 


44.99 


44.86 


0.13 


180 


135 


46.94 


44.93 


0.22 


180 


180 


48-57 


45.00 


0.28 



Table I: The modified mixing angles for various values of phases. Input values are A31 = 
0.002el^^ A21 = O.OOOOSeF^ , = 023 = 45°, ^13 = 0°. 
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ol 


a2 


^12 


^23 


^13 








38.51 


45.00 


0.28 





45 


36.82 


44.93 


0.23 





90 


34.99 


44.86 


0.18 





135 


36.88 


44.93 


0.23 





180 


38.51 


45.00 


0.28 


45 





36.63 


45.06 


0.21 


45 


45 


34.97 


45.00 


0.19 


45 


90 


33.26 


44.93 


0.13 


45 


135 


35.09 


45.00 


0.15 


45 


180 


36.63 


45.06 


0.21 


90 





35.00 


45.13 


0.09 


90 


45 


33.43 


45.06 


0.06 


90 


90 


31.77 


45.00 


0.00005 


90 


135 


33.49 


45.06 


0.05 


90 


180 


35.00 


45.13 


0.09 


135 





36.63 


45.06 


0.21 


135 


45 


35.04 


45.00 


0.14 


135 


90 


33.26 


44.93 


0.13 


135 


135 


35.02 


45.00 


0.19 


135 


180 


36.63 


45.06 


0.21 


180 





38.51 


45.00 


0.28 


180 


45 


36.82 


44.94 


0.24 


180 


90 


34.99 


44.86 


0.18 


180 


135 


36.88 


44.94 


0.24 


180 


180 


38-51 


45.00 


0.28 



Table II: The modified mixing angles for various values of phases. Input values are IA31 
0.002el^^ A21 = O.OOOOSeF^ ^ = 350^ ^ 450^ g^^ ^ qo 
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IV. CONCLUSIONS 



It is expected that symmetries at GUT scale will determine the neutrino mixings. Bi- 
maximal mixing is one of the attractive theoretical ideas proposed in this context. However, 
the solar neutrino data show that the mixing angle 6^12 is substantially less than 7r/4. In this 
paper, we studied how neutrino mass terms arising from Planck scale effects modify the bi- 
maximal mixing scenario. We consider these additional mass terms to be perturbation to the 
main neutrino mass matrix coming from GUT scale and computed the first order corrections 
to neutrino masses and mixings. The changes in all three mixing angles are proportional 
to the neutrino mass eigenvalues. To maximise the change we assumed degenerate neutrino 
masses ~ 2.0 eV. For degenerate neutrino masses, the changes in ^13 and 6*23 are inversely 
proportional to A31 and A32 respectively, whereas the change in 6^12 is inversely proportional 
to A21. Since A31 = A32 ^ A21, the change in 6*12 is much larger than the changes in 6*13 and 
6*23. The actual magnitude of change also depends on the Majorana phases al and a2 but 
is independent of Dirac phase 6, because we assumed that 613 at 0*'' order is zero. In ref. 
[^, it was shown that the change in 613, due to this perturbation, is small. Here we show 
that the change in 6*23 also is small (less than 0.3°) but the change in Ou can be substantial 
(about ±3°). If 6*12 at 0*^ order is 45°, the change induced by the Planck scale terms can 
bring it down to 41.5°, which is the 3a upper limit of the solar mixing angle. 
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